Measuring non-Gaussian fluctuations through incoherent Cooper pair current 
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We study a Josephson junction (J J) in the regime of incoherent Cooper pair tunnehng, capacitively 
coupled to a nonequilibrium noise source. The current-voltage (I-V) characteristics of the JJ are 
sensitive to the excess voltage fluctuations in the source, and can thus be used for wide-band noise 
detection. Under weak driving, the odd part of the I-V can be related to the second cumulant of 
noise, whereas the even part is due to the third cumulant. After calibration, one can measure the 
Fano factors for the noise source, and get information about the frequency dependence of the noise. 

PACS numbers: 74.40.-|-k,05.40.Ca,72.70.-|~m,74.50.+r 



The current in electric circuits fluctuates in time, even 
when driven with a constant voltage. At equilibrium 
or in large conductors, this current noise can be quan- 
tified using the fluctuation-dissipation theorem (FDT), 
which relates the magnitude of the fluctuations to the 
temperature T and the impedance of the circuit. More- 
over, in large wires the current statistics is described by a 
Gaussian probability density which has only two nonzero 
cumulants, the average current and noise power. This 
situation changes for small, mesoscopic-scale resistors ex- 
hibiting shot noise The noise power at low frequen- 
cies is proportional to the average current. Further, the 
statistics of the transmitted charge is no longer Gaussian: 
higher cumulants are finite, and the probability density 
is "skew", i.e., odd cumulants do not vanish. 

For small samples, the frequency scale for the shot 
noise is given by the voltage, eV/H. Shot noise has been 
measured at low frequencies in many types of mcsoscopic 
structures (see the references in However, there 

are only a few direct measurements of shot noise at high 
frequencies uj ^ eV/h and only one of the higher 
(than second) cumulants |j| (at u; ^ eV/h). One of the 
main reasons for the shortage of such measurements is 
the difficulty to couple the fluctuations to the detector 
at high frequencies, or to devise wide-band detection, as 
required for the third and higher cunralants [j|,||5|]. 

In this Letter, we analyze an on-chip detector of volt- 
age fluctuations, based on capacitively coupling a noise 
source to a small J J in Coulomb blockade 0. There, 
the current can flow only if the environmental fluctua- 
tions provide the necessary energy to cross the block- 
ade. In this way, the current through the small J J pro- 
vides detailed information of the voltage fluctuations in 
the source over a wide bandwidth. This information in- 
cludes effects of a non-Gaussian ( "skew" ) environment on 
a quantum system. For the measurement of the second 
cumulant, its characteristics compare well with the other 
suggested on-chip detectors [3,H,S0i0|i based on var- 



ious mesoscopic devices and techniques. The detectors 
proposed in 10, 11] detect the non-Gaussian character of 
the noise, but mapping the output back to the different 
cumulants has not been carried out. The on-chip scheme 
presented here is the first to directly measure the third 
cumulant of fluctuations. In the Gaussian regime, our 
analysis of the noise detection resembles that of Ref. j^] , 
but probes the noise at low measurement voltages. 




FIG. 1: Noise measurement circuit. The voltage noise at point 
A, driven with Vn, induces voltage noise at point B through 
the capacitor Cm and in this way affects J J. J J is in a highly 
resistive environment, R > Rq = /i/(4e^). F^''' are the Fano 
factors for the n'th cumulants in the noise source. The noise 
is read by investigating Im{Vm), as explained in the text. 



We investigate the system depicted in Fig.^ The part 
indicated by the dashed lines represents the noise source, 
and the other part is the detector. The excess voltage 
noise at point A, induced by driving the source (Vat ^ 0), 
gives rise to voltage noise at point B, and the latter can 
be read by examining the current /m(Kn)- There is also 
another source of noise: the equilibrium fluctuations in 
the whole circuit, described by FDT. As shown below, 
the two types of fluctuations can be treated separately. 

In the regime of incoherent Cooper-pair tunneling, the 
I-V characteristics of the detector can be described by 
a perturbation theory in the Josephson coupling energy 
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Ej. This yields for the current through J J [l^ll3| 



neE' 



Im{V^) = [P^(t)(2eK^) - F_^(t)(-2eK„)] . (1) 



Here P^^t)iE) 



J dte'^*^^ ^e">(*)e-"^("') /(27rft) de- 
scribes the Cooper-pair tunneling through the JJ due to 
the fluctuations of the phase difference across the 
junction. This result makes no assumptions about the 
form of the phase fluctuations - only that they do not 
modify Ej itself. As these fluctuations are connected 
to the fluctuations 6VB{t) of the voltage over the J J via 
the Josephson relation, (f){t) — 0(0) + 2e dt'6VB{t')/h, 
Im{Vrn) provides information on these fluctuations. Note 
that, although the coupling capacitance acts as a high- 
pass filter for voltage noise, the conversion to phase noise 
allows to transmit noise down to low frequencies. This 
way the device can be operated non-hysteretically 
and in Coulomb blockade. 

One can identify P,p{E) as the Fourier transform of 
the moment-generating function k^{x) = (e-^'^'^'^e"-^'^'-''-') 
of — (/)(0), evaluated at x = *• This can be expanded 
in the cumulants C^(t), ^^(x) = exp[xV(2!)C|(t) + 
X'^/(3!)C3 (t) -|- •••]. These cumulants are defined such 
that 0(t) is ordered before (j){0) in the expectation value. 
The expansion defines a function J0(<) = ln(K0(i)) = 
J2 (t) + Jsit) + ■ ■ ■ , where for stationary fluctuations 

J2(i) = mt) - mm)) and Mt) = ^mmt) - 

(j){0)](f>{0)) , etc. In the Gaussian hmit, J^ft) = J-,j,{t) = 
•hit) coincides with that applied in Ref. 

The odd cumulants of the fluctuations in the source 
break the symmetry between the positive and the nega- 
tive Vm- To separate the non-Gaussian effects, we con- 
sider the even and odd parts of the I-V, /s/A(Kn) = 
(^m(Ki) ±-^m(-Kn))/2. The odd part /^(Ki) describes 
the general behavior mostly due to the even cumulants 
(c.f., Eq. © below), and the even part /s(Kn) (Eq. Q) 
responds to the odd cumulants, vanishing if C2n+i = 0. 
We show below that tuning the voltage Vm and measur- 
ing /A/s(Kn) gives access to the frequency dependence 
of the lowest cumulants at large bandwidths. 

We have to take the additional Gaussian fluctuations 
Jjp{t) from the total impedance of the setup into account, 
as they will inevitably influence the measurement. Thus, 
we split the fluctuations as parts J0(<) = J^(t) + J^{t) 
with the excess fluctuations due to the driven source de- 
noted by J^{t). The latter includes the non-Gaussian 

effects for which jf^ ^ J^. We also define P^^^{E) = 

J dte'^^l^{e^"'''^^^)/{2-Kh). Using detailed balance [l^, 
we obtain for the current lD{Vm) = -fm(Kn; = 0) = 
■neE'][Pf{2eVyn){l - cxp(-/32eKn))]/;i. In the presence 
of the excess noise, the current is given by 



1 - p-PE 



„/3(iS-2eV„) pS 



evaluated at -B = 2eVm- Here = ksT and * denotes 
the convolution over the energy. Thus, the detector char- 
acteristics, described by P^{E), do not need to be known 
exactly, but they can be calibrated by measuring the cur- 
rent iDiVm) in the absence of the additional noise. 

For connecting the fluctuations at the JJ and the 
source, we relate the mutually uncorrelated intrinsic fluc- 
tuations Sla, Sib and Sim through the resistors Ra, 
Rb and R to the voltage fluctuations AV at point B, 

AViLu) = R[-SIra{uj) " iLoCraRs{SIb - Sla - SIm)]/G{Lo) 

through circuit analysis. Here G(w) = 1 — iuj{RC + 
RsCm) -u;^RCjRsCm, C = Crn + Cj and Rs = 
{R^^ + Rb^)~^- The equilibrium fluctuations in the 
source are present even for Vjv = 0, and they can be 
included in the calibration through the FDT. The ex- 
cess fluctuations due to driving produce excess phase 
fluctuations, characterized by the n'th order correlators 
((/)(wi) • • • (/)(w„)) = 27r(5(wi H h UJn)Sn^iuj) [11, 



Sn^{uj) = A" 



(3) 



where Co 



V(4e2), 



(wi, . . . ,a;„), A : 
= RC, G„(cD) 



TrRsCm/{RQC), Rq = 
= G{uji) ■ ■ ■ G{LL>n), and 
{SI{ll!i) ■ ■ ■SI{ujn)) = 2tt5{uji + . . . uJn)Snii^)- Thus, we 
find that the excess phase noise and its cumulants are 
governed by powers of A, the current Int/b and the band- 
width of the current correlators, times r. In the absence 
of the source {Cm 0) , the detector phase fluctuations 
Jff/7 (t) have been calculated in (see also Fig. 11 of 
|l2|). The equilibrium fluctuations in the source slightly 
modify this behavior ^3- ^ RsGm ^ RCj, the result- 
ing J^{t) is close to JRc{t), with the capacitance given 
by Gm + Gj and resistance by R, hence the JJ is in the 
insulating state for R > Rq. Below, we concentrate on 
this insulating parameter regime. 

Now assume the source is driven weakly, such that 
J^{t) ^ 1. In this case, one may expand the exponential 
e-'^W « 1 + J^(t) [13. Then, 

P^iE) = 5iE)+^^-S,,iO)SiEy' '''^''"'^ ' 
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Here S2rf>{uj) is the driven phase noise spectrum induced 
by the source, S2(f,{t) its inverse Fourier transform, and 

1 

K^{U)) — — j duj'lTCLS^cl,{u)TU)' ,—U) — Lo') (5) 
I" J-oo 

describes the third cumulant of phase fluctuations. In the 
derivation of this form for K^, we used the hermiticity 
of (j){t) and the stationarity of the fluctuations. In what 
follows, we cut the expansion in the third cumulant |l5l |. 

Let us first concentrate on the antisymmetric part of 
the detector current. In the first order in J^{t), it only 
probes the even cumulants. In this case, it can be ex- 
pressed as lA{Vm) = loiVm) + SlA{Vm), where for sym- 
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metric S2i{uj) = S2i{-uj) [3, using Eq. ®, 

/OO 7 
^Va{u;; V„,){S'M + S^Aco)). (6) 
-OO 

Here 2?a(cc';K„) = \^T^[ID{V„^ + hw/{2e)) + loiVm - 
huj/{2e)) - 2/D(\4t)]/(2e^G'2(w)) characterizes the fre- 
quency band for the detection of 5*2/(1^). This band can 
be tuned by tuning the bias Vm (sec Fig.|21l. 




FIG. 2: (Color online) Left: characteristic function 
T>A{u};Vm.) of the excess noise measurement, (c.f, Eq. ®), 
probing the spectrum S2i{i^) of current fluctuations. Dot- 
ted line shows the frequency dependence in 821(^0) (arbi- 
trary units), for T — and V = e/C. Right: function 
[loiVm + hu!/{2e)) — loiVm — fttj/(2e) )] /2, characterizing the 
measurement of the third cumulant. The frequency depen- 
dence in K^{lj) (arbitrary units) is shown in black. In the 
limit eVn ^ fi/T, the width of K^{u) is given by 1/r. Here, 
Ec = 2e^lC, R = 6i?Q, Ra = Rb = 0.1 Rq, and'c„ = WCj. 
Changes in Cm do not essentially affect the figure. 

The even part of the detector current can then be re- 
lated to the odd cumulants of phase fluctuations. For a 
weakly driven source, assuming K^{'-uj) ~ —K^{lj) (see 
below), it is given by 

Is{V^) = J/z. (v„, - ^) K, (CO) . (7) 

In this way, Is{Vm) probes the frequency dependence of 
the third cumulant of source fluctuations, (c.f., Fig.EJ. 

Let us now analyze 5lA{Vm) for the noise source shown 
in Fig. ^ The frequency dependence of the nonsym- 
metrized noise S'/(w) is derived, e.g., in Ref. @. This 
derivation holds provided that the Thouless energy Et 
for the resistors greatly exceeds cVn- Subtracting the 
fluctuations at Vn = we get the excess noise S2i{'jj) = 
FiS{Vi)/R„ where i = {a, 6}, = R^VN/{Ra + Rb) and 

el^sinh(f)-2fi^coth(^)sinh^(^) 
^ ' cosh(f^)-cosh(|f) ■ ^ ' 

An example of the frequency dependence of iS'2/ is plotted 
in Fig. El Note that it is symmetric with respect to the 
sign reversal of w. This property can be traced to the fact 
that in our example the source impedance, characterizing 
quantum fluctuations, stays constant. Now, S21 can be 
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FIG. 3: (Color online) Antisymmetric 5lA{Vm) (solid, left 
axis) and symmetric part Is{Vm) (dashed, right axis) of the 
detector current variation. Dotted lines are expansion results 
(Eqs. (|7|l). Parameters are as in Fig. |5| The source is 
assumed to consist of a mesoscopic {F^ ^ 0) and a macro- 
scopic resistor (_F,j — 0) for which Sni{<^ = 0) = e"~^_F„Jjv. 
Inset: 5lA{Vm = e/2C) (solid) and Js(V„ = e/2C) (dashed) 
as functions of Im- Scaling parameters are D = F2DA = 
eF2 J dcLiI'A(i^; Vm)/(27r) for the antisymmetric part, and 

D = F-^Ds = Fi^^X" j duIoiV^ - huj/2e)Lj[4 + Sicjrf + 
(tjr)'']^^ for the symmetric part. 

substituted to Eq. ^ to find the effect of the driven 
noise on the detector current. Figure 13 shows an example 
of SlA{Vm), i.e., the probed shot noise, for a few bias 
voltages Vn- For Vn » h/r, shot noise is essentially 
white over the detector bandwidth (Fig. [SJ and thus the 
signal is linear in V^v- In this case, lA{Vm)/lN (Fig. |3| 
inset) depends only on the factor DAiYm), which can 
be related to the calibration measurement, and the Fano 
factors which can thus be measured from this curve. 

Next, consider the detection of the third cumulant 
of current fluctuations. Its frequency dependence has 
been described in various limits for different systems in 
Ref. At T = 0, the third cumulant at zero frequency 
is of the form ^3/ = F^c^In, and its dispersion occurs on 
the scale lUc = min^eV/h, Ex/fi, I/e). For RsCm ^ RCj 
and a frequency independent S31 within the detection 
bandwidth, the resulting K^(lu) at T = is given by 

(9) 

This is an antisymmetric function of w, due to the simple 
form assumed for 5*3/, and its frequency scale is given 
by 1/t. With the knowledge of the detector calibration 
current /_d(14i), plugging K^{uj) into Eq. iQ allows us to 
calculate the response of the symmetric detector current 
-^5(^71) to the third cumulant in principle for any type 
of resistors Ra and Rb- A few examples of Is{Vm) are 
shown in Fig. |3 As for the second cumulant, following 
Is{Vm)/F)s{Vm) allows to measure the Fano factors. 

The general measurement scheme is valid at a finite 
temperature as long as the Coulomb-blockade condition 
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FIG. 4: (Color online) Antisymmetric (solid) and symmetric 
(dashed) detection currents 5JA(Kn) and /s(Kn) for the same 
parameters as in Fig.|31but with a few different temperatures 
at fixed Viv = "i^/C, i.e., Im ~ 60e/r. For Is{Vm), we neglect 
the temperature dependence of the intrinsic fluctuations. 



T < Ec/kB is satisfied. The effects of T > are il- 
lustrated in Fig. ^ Compared to the previous example, 
two corrections arise: the function /^((Kn) characteriz- 
ing the detector becomes smoother and its amplitude de- 
creases, and the form of the source excess fluctuations 
changes. The two are characterized by different temper- 
ature scales, Ec/kB and eV^ /kB, respectively. 

Measuring the calibration current loiYm) along with 
the antisymmetric and symmetric currents, SlAiYm) and 
IsiVm), allows to find the second and third cumulants of 
excess current fluctuations in the source within the band- 
width described in Fig.|21(for typical parameters 
the range of 100 GHz) to an accuracy limited mostly 
only by the resolution of the current measurement (res- 
olution of 0.1 pA yields ~ 100 (fA)^s for the second and 
- 0.01 (fA)3s2 for the third cumulant 0, [13). In the 
limit Cm ^ Cj, RsCm <C RCj, the only information 
required about the setup are the resistances i?Q, Rb and 
i?, and the sum capacitance Cm + C'j, all of which can 
be measured separately. Thus, the scheme allows an ac- 
curate determination of the second and third cumulants 
over a large bandwidth and for the third cumulant, over- 
comes the bandwidth problems encountered in Ref. 0. 

In summary, we show in this Letter how incoher- 
ent Cooper-pair tunneling in small JJs with a high- 
impedance environment can be used for accurate and 
wide-band detection of voltage fluctuations. Via the sym- 
metry of the detector output current, one can identify the 
contributions from the second and third cumulants sep- 
arately. While the presented example is on the measure- 
ment of fluctuations in samples exhibiting no Coulomb 
blockade, this is not a limitation of the scheme itself. 

Note added in proof: During the refereeing process, we 
became aware of Ref. , which points out that the Fano 
factors for the third cumulant depend on the deflnition 
of the measured observable. Given that F^ are related to 
the nonsymmetrized observable, our results remain valid. 

We acknowledge the useful discussions with W. Belzig, 
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